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Abstract 

We identify a fractal scale s in a family of Borel probability mea- 
sures ^ on the unit interval which arises idependently in quantum 
information theory and in wavelet analysis. The scales s we find sat- 
isfy s £ M+ and s ^ 1, some s < 1 and some s > 1. We identify these 
scales s by considering the asymptotic properties of u (J) / \J\^ where 
J are dyadic subintervals, and \J\ — > 0. 

I Introduction 

While finite non-commutative operator systems lie at the foundation of quan- 
tum measurement, they are also tools for understanding geometric iterations 
as used in the theory of iterated function systems (IFSs) and in wavelet anal- 
ysis. Key is a certain splitting of the total Hilbert space and its recursive 
iterations to further iterated subdivisions. This paper explores some impli- 
cations for associated probability measures (in the classical sense of measure 
theory), specifically their fractal components. 

In quantum communication (the study of (quantum) error-correction 
codes) , certain algebras of operators and completely positive mappings form 
the starting point; see especially the papers [18] and [21]. They take the 
form of a finite number of channels of Hilbert space operators Fi which are 
assumed to satisfy certain compatibility conditions. The essential one is that 
the operators from a partition of unity, or rather a partition of the identity 
operator I in the chosen Hilbert space. Here (Definition II.1|) such a systems 
[Fi) are known as column isometries. An extreme case of this is when a 
certain Cuntz relation (Definition II. 1|) is satisfied by. Referring back to our 
IFS application, the extreme case of the operator relations turn out to cor- 
respond to the limiting case of non-overlap. Using this operator theory, in 
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this paper we explore the fractal measures associated with the inherent self 
similarity affine fractals, a subject involving both iterated function systems 
(IFSs), and an aspect of quantum communication. 

In the paper we aim to draw up connections between the following three 
areas, quantum channels, fractal measures, and wavelets. Our understanding 
of quantum measurement follows the tradition of Kraus [20]: By a quantum 
measurement of a system refers to a system in some state and we want to de- 
termine whether it has some property E, where E should be thought of as an 
element in a logic lattice of quantum yes-no questions. Measurement means 
submitting the system to some procedure to determine whether the state sat- 
isfies the property. The reference to system state in turn must be given an 
operational meaning by reference to a statistical ensemble of systems. Each 
measurement yields some definite value or 1. 

The connections to fractal measures and wavelets (see e.g., [16]) derive 
from a common mathematical core based in turn on hierarchical patterns 
common to them. Wavelets are computational bases in Hilbert space which 
make use of the scaling law, with each step in the scaling creating a refinement 
in a resolution of data, and taking advantage of similarity from one level 
of resolution to the next. Hence we arrive at the characteristic feature of 
fractals, see [1], [2]. 

Mathematically, an identification of quantum channels may be made with 
a finite system of operators (Fj) which together form a column isomety; see 
(1) below. In the physics literature, they are known as Kraus-systems, and 
mathematically they generalize the so called Cuntz relations; see |2], 
[II, [23], [5], [12], [20]. The systems (F,) which are the focus of our paper 
in turn determine completely positive mappings a as recalled below. In 
this paper we show that repeated application of operator system yields the 
kind of refinement that underlies both wavelets and fractals, and moreover 
that the statistical properties from quantum channels is of significance in 
decomposition problems from wavelets. And with the scaling from wavelet 
theory now corresponding to a. 

Both the quantum mechanical measurement problem and IFSs have as 
starting point a finite set of operations: in the case of IFSs they are geo- 
metric, and in the quantum case, they involve channels of Hilbert spaces 
and associated operator systems. The particular aspects of IFSs we have in 
mind are studied in [10]; and the relevant results from quantum communi- 
cation in [18] and [21]. We begin the Introduction with some background 
and motivation on IFSs. The operator theory, the fractal measures and their 
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applications are then taken up more systematically in section II below. 

Let 7i be complex Hilbert space, and let A be a finite set. We will be 
interested in an indexed set of operators {Fi\i G A} satisfying 

where / denotes the identity operator in Ti. 

Definition I.l ^ finite system of operators Fi in a Hilbert space Ti is said 
to be a column isometry if 



e n 



is isometric; and it is said to be a Cuntz-system if also 

FiF* = SijL 

Because of a certain reasoning outlined in the references below such sys- 
tems are called measurements in quantum probability; see e.g., [16], [2H 
OA0404553], or [21], [25l quant-ph 101061], [S].But they arise in other fields 
as well, in representation theory, in geometric measure theory, and in wavelet 
analysis; see e.g., [8], [H], [13], [15], and [HI Minimality, Adv. Math]. Closely 
related systems of Hilbert space operators play a big role in the theory of 
frames [TT], [3] and their engineering applications. Suppose = N. De- 
note by N the cyclic group of order viz., Z/NZ = {0, 1, ■ ■ ■ , A^ — 1} or 
Ztv, and set 

= = {all functions N ^ Zjv} . (2) 

We shall give Q its Tychonoff topology, and we view it as a compact Hausdorff 
space. 

Let JF = (Fj)jgj^ be a measurement, and let G 7i be a unit- vector, i.e., 
a quantum mechanical pure state. Then it is immediate that 

fi^{t): = \\Fiij\\\ zGN (3) 

is a probability distribution on N. 

The measure in the title of the present paper refers operator-valued, or 
scalar- valued measures on Q = induced by ([1]) - ([3]). the induction from 
Ztv to Q]\r is based on the Kolmogorov consistency condition [19], [H], as 
follows: 
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Definition 1.2 

(a) Cylinder sets: Let (ii, ■ ■ ■ , ir) £ and set 

C (zi, ■■■ ,ik) = {w e Vt\w (1) = zi, ■ ■ ■ , w (/c) = ifc} . (4) 

(b) Operator valued conditional probabilities: 

P{C{i,,--- ,i,))=FI...F:F,,...F,, (5) 

(c) Kolmogorov consistency: Formula below. 
Since we have the disjoint union 

C {ii,- ■ ■ ,ik) = IJ C(ii,--- ,Zfc,j); (6) 
we need the formula 

P(C(2i,---,^,))= J]P(C(2i,---,^fc,j)) (7) 

in order to extend ^ to a probability measure P defined on the Borel subsets 
ofQ. On the other hand, it is easy to see that © is satisfied by Just 
use the basic formula for the given measurement T . 

Lemma 1.3 There is a unique positive operator-valued probability measure 
P defined on Vt, and satisfying @ For any Borel set B d Vt, P (B) is well 
defined, and {ip\P (B) ip) > for all ip E 7i. Moreover, P is sigma- additive, 
or countably additive, i.e., 

oo 

whenever Bi,B2, - ■ ■ are disjoint Borel sets; 

p{n) = i (9) 

Proof. The argument for the existence and uniqueness of the extension is a 
standard application of Kolmogorov consistency. See \X9l, ^14j, ^15j, \X3\ for 
more details. ■ 
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Examples 1.4 

(a) N = 2, n = C^, ^= (J), and 

Then 

fi^{C{t,,--- ,ik)) = 2-'', I, e {0,1}. (11) 



(b) N = 3, H = C^, ip 








Then, for ij G {0, 1,2}, G N, and 1 < j < fc, we h 



ave 



H (C* (^1, ■ ■ ■ ,4)) 



if some ij = 
otherwise. 



:i2) 



(13) 



Proposition 1.5 If we introduce N-adic partitions of the unit-interval [0, 1] 
as follows 



C{^^, 



.Ik) 



ii ik ii ik 1 
_!. J \ ^ + \ ^ _] 



(14) 



then the measures P {■) and fi^ 



|P (■)'?/') are Borel measures, each 



supported on [0,1]. If N = 2, the measure p7|) in Example L4_ (a) turns 
out to be merely Lebesgue measure restricted to [0, 1]. If N = 3, the measure 
(flgj) in Example \I.4\ (b) is the middle-third Cantor measure supported in the 
Cantor set X3. 

Remark 1.6 Recall ^3 is the unique (compact) subset o/M satisfying 

3X3 = X3 U (X3 + 2) , (15) 
and the Cantor measure fi = fi^ is the unique Borel measures satisfying 

'x + 2' 



I fix) d^{x) = \(^j f (I) d^{x) + j f 
for all bounded Borel functions f . 



djj, {x] 



(16) 
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Proof, (of Proposition 1.4) The assertions follow from standard appli- 
cations of the Kolmogorov extension principle, and the reader is referred to 
[15] . [13] for additional discussion. ■ 

The column isometrics (Fi) introduced above can be viewed as Kraus 
operators [20] from the theory of quantum channels, modeling "instruments" 
in Kraus's formulation, and operating on quantum systems by producing a 
classical measurement (outcome.) The theory of finitely correlated states 
[12] has an instrument generating a classical state which takes the form of 
a measure on an infinite product of a finite alphabet; a construction which 
parallels the theme of our paper. It would be intriguing to explore how the 
physics of [12] refiects itself in the measures more directly associated to fractal 
measures, and wavelets. Part of the answer lies in how families of wavelet 
packets adapt to signals (or in 2D) to images, see e.g., [16]. The operators 
in the system (Fj) are iterated in steps with each iteration step creating a 
subdivision of the masks in the previous more coarse resolution. A monomial 
of degree k in the generators Fi corresponds to a fold subdivision. Even 
though this algorithm is "classical," the generators Fi are non-commuting 
operators. The limit of this iterative scheme as k tends to infinity is made 
precise by the measures that we analyze in theorems II. 3 and III. 2 below. 



II Fractal Scales 

The authors of [9J recently adapted the discrete wavelet algorithms to fractals, 
and the present work extends [9]. 

The distinction between the two prototypical cases (a) and (b) in Example 
11.41 can be made precise in a number of different ways; for example, it can be 
checked that the fractal dimension (in this case = the Hausdorff dimension) 
of (a) is 1, and of (b) it is s = |^ = log3(2). For our present discussion, the 
following definition of the fractal dimension will suffice: If a subset X C M"^ 
is obtained by the iteration of a finite family of contractive and affine maps 
M — i> M then the fractal dimension s of X 

log {number of replicas) 
log {magnification factor) 

Following Proposition 11.51 especially ( |T4|) , we will restrict attention in the 
following to subsets of [0, 1] and measures defined on the Borel subsets of 
[0, 1]. If J C [0, 1] is a subinterval, we denote by \J\ the length of J. 
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Definition II. 1 Let fi be a probability measure on [0, 1] defined on the Borel 
sets. 



We say that s_ is a lower scale of if 



liminf(^ > 0; (18) 

and we say tliat s+ is an upper scale of if 

iimsup . J.S+ < oo. (19) 

It is easy to see that the Cantor measure = /i^ in Example [L4] (b) has both 
upper and lower scale s+ = s_ = s = logs (2)- 

Our next result is motivated by examples from wavelet analysis. Before 
stating our general result we first recall the wavelet examples. To emphasize 
our point, we do not consider the wavelet examples in the widest generality. 

Example II. 2 Discrete Wavelet Transforms. Let ao,ai, - - be a se- 
quence of complex numbers such that 

'^a]aj+2k = So^k- (20) 
j 

In the summation ( \2(M . it is understood that terms are zero if the subindex 
is not in the range where a ^ 0. 

We define operators Fq and Fi on the Hilbert space T-C=i'^ as follows: 

(^oOr = 5Z«2,-fca (21) 

k 

and 

k 

Then it is easy to check that ([1]) holds, and so the pair [Fq, Fi) defines a 
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measurement in the sense of the definition in Section I. In this case, more is 
true: The adjoint operators F* are isometries with orthogonal ranges, i.e., 

FiF* = S,,,L (23) 

If the sequence oq, Oi, ■ ■ ■ from (I2U]) is finite, then it is easy to see that the 
number of non-zero terms must necessarily be even. We consider 2D scalars. 



ao, ai. 



(24) 



and the corresponding two {2D — 1) by {2D — 1) matrices Fq and Fi defined 
as follows: 



a2 





ai 
as 




ao 
a2 



\ 


















ai 

02 



(12D-2 












ai 

03 



\ 







ao 

a2 



•• 

and Fi built the same way, but using the numbers 

bk- = {-I)'' a2D-i-k 
For D = 2, the two matrices are simply 



ao 
-fo = I ai tti I and Fi 

a3 a2 



(25) 



(120-1 (120-2 J 



(26) 















ai 







Oi 



(27) 



Staying with ao,ai,a2,a3, there are practical reasons in wavelet analysis to 
add the following two requirements to ( !20l) : 



1 ai e 



and 



:ii) J2a, = V2. 

i=0 



Taking the combined conditions together, it can be shown that 00,01,02,03 
are determined by a single real parameter (3 thus; 



ao = ^ (1 + y2cos/3) 
ai = ^ (1 + 72 sin/?) 
«2 = ^(1- V2COS/5) 
^3 = ^(1- v^sin/3) . 



(28) 



A consequence of (128|) is that each of the three pairs (oo,oi), (00,03), and 
(oi, 02) lies on the circle 



X 



1 V 



1 



2^2 



1 

4' 



(29) 



see Fig. 1. 




Figure 1. One of the three pairs (oq, 03) 



Theorem II. 3 Let the numbers oq, o„, ■ ■ ■ , 02D-1 be given, and suppose 

is satisfied. Let Fq and Fi be the corresponding matrices determined by ( f^3|) - 

^2Bi . Suppose further that oq ■ 02_d-i 7^ 0. Let 

a: = max (|oo|^ , |o2D-i|^) ■ (30) 
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Then the number 



is a lower scale of /iq 



In a 



In 2 



(31) 



(col-P (■) Co) where 

( 1 \ 





eo 



(32) 



It follows in particular that if a > ^, then /io has a lower scale s > 1. 
Moreover, we can get the lower estimate i\ lS\) satisfied for dyadic intervals 
J inside any non-empty open subset of [0, 1]. 



(33) 



Proof. From the representation (!25l) of the two matrices Fq and Fi we 
conclude that 

F*eo = a^eo and F*eo = 020-1^0. 
With the dyadic representation f |T^ . N = 2, set 



and 



_L J L _1 

2 2^= 



F. . . . F. 



(34) 
(35) 



Then it follows from (fTSll and Lemma 11.31 that 

/^o + 2-'')) = {eo\F*^ ■ ■ ■ F:F,^ . . . F,,eo) 
= {Fi^ ■ ■ ■ Fi^eo\Fi^ ■ ■ ■ Fi^eo) 
= WF^eof 
>\{eo\F^eo)f 

(Schwarz) 

_ I |2.#(i=0) I |2-#(i=l) 
— I^ol ■ \0.2D-l\ 

using 

Let \^ be a non-empty open subset of [0, 1], and pick k G Nand as in 
such that the interval ^ + 2~^^ is contained in V. 
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We now turn to the two possibilities for the number a in fl30p . 



2 

If a = |ao| , then 

/io([e,e + 2-'-")) = ||Fo"F5eo||' 

and we conclude that the expression 

/io([e,e + 2-^-")) 
2" s(fc+?i) 



(36) 



is bounded below as — > oo, and hence f|T8l) holds for s = log2 (a see 



2 

If instead a = \a2D-i\ , then 



2#(i=0) ^#(i=l)+. 



/xo([e + 2-'=(l-2)-"),e + 2-'=) > |ao| a 

by the same reasoning used in the first case. We now get the lower estimate 
(fT8|) satisfied for the intervals J = + 2"'= (1 - 2^") , ^ + 2"'=) as ^ oo. 
This completes the proof ■ 



III Upper and Lower Fractal Scales for the 
Measure /xq 



Consider the example outlined in (1281) above. The two matrices Fq and Fi 
are used in wavelet analysis where they refer to low-pass and high-pass filters; 
terms that derive from signal processing, see [6] and [3]. 

Recall that when 00,01,02,03 are given by fl28l) then there are solutions 

(j),^p in (M) to 

(f){x) = V2 (oo0 (2x) + oi0 (2x - 1) + O20 {2x - 2) + 03^ (2x - 3)) 
Jjg (x) dx = 1 

ij{x) = V2 (o30 (2x) - O20 {2x - 1) + oi0 (2a; - 2) - Oo0 (2x - 3)) 
/m ^ (a;) rfx = 0, 

and when /3 G ]R\ ({±^,±^} + Z27r) , then the two functions (the scaling 
function) and (the wavelet) satisfy the further conditions 
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/ (j) (x) (f){x — k) dx = (5o,fc5 k E Z 

and 

{2^/2^ {2^x - k) \j,k eZ} 
is an orthonormal basis for the Hilbert space (M ) . 

The general feature of the matrices ( 1271) is the slanted shape; and the FiS 
are called slanted Toeplitz matrices. 

Example III.l A one-parameter family of wavelets. For the case D = 
2 treated in ( Tggj) . we have 

'''' 

and it follows that the number 

a = max (oq, 03) (38) 

satisfies a > | when the parameter P is in one of the two intervals: 
(0 \P\ J where a = a^; 

or 

(a) — ^ < P < —J where a = al; 
see Fig. 2 below. In these two regions the scale number s satisfies s > 1; 
hence the fractal feature of the measure hq. 
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Figure 2. The two functions = Oq and 03 = as fl28l) . Two 
regions (i) and (ii) with scale number s satisfying s > 1. 



When the four numbers oq, cii, ^3 that make up the matrix Fq are given 
by formulas (128|) . one easily computes the spectrum of Fq as follows: 

/ ^ sin/3-cos/?\ 
spec(Fo) = |ao(/5),--^, ^ j>, (39) 

and we have sketched the point 

A (/?) = 'M^^ (40) 

in the spectrum in Figure 3. An inspection shows that A (/3) is not dominant 
in spec (Fq) in the sense that the inequalities 

ao(/3)>-^>|A(/3)| (41) 
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hold in the region (i) from Figure 2. 



1+^ 




4^) 




The eigenvalue A = 2Mz£2^ jg not dominant. (i):ao > -^^ > A {pf 



We will now turn to our analysis of the two-sided scale bound for the 
measure /xq aiid we show how it applies to the matrices Fq and Fi which are 
used in wavelet theory. 

Theorem III. 2 Let the numbers ao,ai,--- ,a2_D-i satisfy condition {\2U[\ . 
and in addition 



j=0 

The two matrices Fq and Fi are defined as in - (f^6j) . We make the 
following additional assumptions on the spectrum of Fq: 

{i) ao ■ a2D-i 7^ 0; 




(42) 



(43) 



(a) \ao\ > max {|A| A G spec (Fq) \ {ao}} ; 



(44) 



(Hi) the algebraic multiplicity of ao mspec(Fo) is one. 



(45) 
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Then there is a unique vector v such that 

Fqv = aov and {eo\v) 



(46) 



where Cq is the vector i\3^ with 1 is the first slot and zeros in the rest. 
Moreover 

. = log, (W-) = -il^ (47) 

is both an upper scale and a lower scale in every non-empty open subset of 
[0, 1] for the measure 

/io(-) = (eo|P(-)eo). (48) 
Proof. Set w = ( 1, 1, - ■ ,1 ). Then implies 

2D-1 ones 

1 1 

wFq = —/=w or equivalently FqW* = —/=w* (49) 
V 2 V 2 

where w* denotes the column vector corresponding to w. So property (jS]), 
i.e., (ii) above, yields 

|ao|^ > - and |ao| > |a2D-i| • (50) 
2 

Since -FQ*eo = aoCo, ( H5ll implies that is not in the spectrum of the matrix 
G arising from Fq by deletion of the first row and the first column, i.e.. 



as 



G 



ao 

02 





ai 





ao 











0-20-1 02D-2 0-20-1 O2O-I 




V 





as a2 



a2D-i 020-2 / 



(51) 
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Now define a vector v in ^ by 



W = Co + (ao/2D-2 - G) 



-1 



/ "2 \ 

a4 







Vo / 

To better visualize (ISTIl the reader may check that, if D = 2, 



ao 

^2 Ctl CtO 

03 a-2 



G 



(X\ (2q 
03 02 



and 



(1\ CLq 

a-3, a2 



a2 




/ 1 

p{ao) 
03^2 
\ p(ao) 



where 



p (ao) = cto ^ (cti + 02) «o + ^102 - aotts. 
Recall that the characteristic polynomial of G is 

p (A) = - (trace G)X + det G, 



(52) 



(53) 



(54) 



(55) 



and that p (ao) in the fraction of is evaluation of at A = ao- Hence, 
assumption comes into play. 

Returning to the general case, we claim that v satisfies ( H6|) . Indeed, let 
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V be given by fl52p . then 



Fv = aoeo + {l + G {a^I - G)'^) 



04 







aoeo + clq {a^I - G) ^ 



Vo 

04 



0-20-2 





= flow, 
and 

{eo\v) = 1. 

This proves fj46|) . Conversely, using fH^ - fj45|) . one shows that the solution u 
to fH6l) is unique. 

Now let y be a non-empty open subset of [0, 1]. Then pick k G Nand a 
dyadic fraction 



^ = :± + . . . + (56) 

such that i + 2-'=) C V". 

We now wish to estimate /iq ( ^ + 2"*^"")) and get the asymptotic scal- 
ing rate as n ^ 00. 

To that end, we prove in the next section (in a separate lemma; see espe- 



cially ([59]) that lim^^oo ^^'FoF^eo = {eo\F^eo)v = {F^eo\eo)v 



a, 







^2D-1 



and result lim |ao| "/io([^,^ + 2 ")) 
cxD. Since ||f ||^ > 1, the desired conclusion follows. ■ 

In applications to wavelets, the measures in the title of the paper are 
used in the computation of transition matrices for transformation between 
two orthogonal families in the Hilbert space (M) : 



2.#(i=0) I |2.#(i=l) II ||2 

Oo |a2D-i| \\v\\ n 



17 



(i) a wavelet basis (2^/^^ {2Px - k))^^^^, 

and 

(ii) a wavelet packet (0„)„gN^^ Nq = {0, 1, 2, ■ ■ ■ }; 

02n (a;) = V^Y.k& 0n (2x - k) , 

02n+i (a;) = V2 Zlfcez ^fc ^ri (2x - A;). 
The adjustment of dyadic scaling in (i) is made with variations in p G Z; 
and hence with the size of the dyadic intervals J {k,p) = [k2~'P, {k + 1) 2^^). 
The concentration of mass at each J {k,p) is determined by the measure. 



IV A Technical Lemma 

In the proof of Theorem IIII. 21 above, we relied on the following lemma regard- 
ing operators in a finite-dimensional Hilbert space. While it is analogous to 
the classical Perren-Frobenius theorem, our present result makes no mention 
of positivity. In fact, our matrix entries will typically be complex. 

Notation IV. 1 If Ai is a complex Hilbert space, we denote by L{Ai) the 
algebra of all bounded linear operators onM.. IfA4 is also finite- dimensional, 
we will pick suitable matrix representations for operators F:A4 A4. 

Suppose J\4 contains two subspaces A^j, i = 1,2 such that Aii ± M.2 
and Ai = Ail © A42, then we get a block-matrix representation 

where the entries are linear operators specified as follows. 

A:Mi Ml, B:M2 ^ Aii; 

and 

C:Mi M2, D:M2^M2. 

If dim All = 1, and Aii = Cw for some w & At, then we will identify the 
operators Ati At with Ai via T^:C 9— >■ zt], where f] G Ai. the adjoint 
operator is T*x = {rj\x)w, x G Ai. 
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Lemma IV. 2 Let M. he a finite- dimensional complex Hilbert space, with 
d = dimA^. Let F G L{Ai), and let a ^ C satisfy the following four 
conditions: 

{i) a G spec (F) ; 

[a) \a\ > max {|A| |A G spec (F) \ {a}} ; 

(Hi) the algebraic multiplicity of a is one; 

{iv) there is a w E M., \\w\\ = 1, such that F*w = dw. 

Then there is a unique ^ G 7W such that 

= 1 andF^ = a^. (58) 

Moreover, 

lim a~'^F'^x = {w\x)^ for all x e M. (59) 

n— »oo 

Remark IV. 3 There is a constant C independent of d = dimA^ and x, 

such that 



{ 


s 


n 




a 





spec(F)\{a}| (60) 

Proof. (Lemma iv;!) Set M^:= MqCw = {x E M\{w\x) = 0}. Then 

M = Cw®M, (61) 
and we get the resulting block-matrix representation of F, 

a 00 ■■■0 



T] G 



(62) 



where a is the number in (i), the vector r] G Ai'^, and operator G E L (A^"*") , 
are uniquely determined. 

As a result, we get the factorization 

det (A - F) = (A - a) det (A - G) (63) 

for the characteristic polynomial. Assumptions (ii) and (iii) imply 

spec (F) \ {a} = spec (G) ; (64) 

and in particular, we note that a is not in the spectrum of G. Hence the 
inverse (a — G)-^ is well defined, and (a - G)~^ E L [M^]. 
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We claim that the vector 



= w + {a — G) ^ 1] 



(65) 



satisfies the conditions in (1^ . 

First note that {a — G)^^ rj G ^A^^, so {w\C,) = {w\w) = \\w\f = 1. 
Moreover, 

= aw + T] + G {a - G)'^ r] 
= aw + a{a — G)^^ rj 

which proves the second condition in fISH]) . Uniqueness of the vector ^ in fl58l) 
follows from flMl) . 

Using the matrix representation (l62l) . we get 



/ 2 



a?] + Grj 



00 ■••0 \ 
00 • ■ ■ 



QT 



J 



and by induction, 



a" 00 ■ ■ ■ 

^n-i^ + a^-^Gr/ + ■ ■ ■ + G^'-^T] G" 

a„ 00 ■ ■ ■ 

{a"" - G'') {a - G)~^ rj G" 



(66) 



Hence, if we show that 



lim a-"G" = 0, 



(67) 



then the desired conclusion fl59l) will follow. Using the matrix form (l66l) . the 
conclusion (|59|) reads 



lim a-"F' 

n— >oo 



(a-G)"S 



00- --o 





(68) 



In proving fl67jl . we will make use of the Jordan- form representation for 
G. Jordan's theorem applied to G yields three operators D, V , N E L (A^"*") 
with the following properties: 
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(1) D is a diagonal matrix with the numbers spec (F) \ {a} down the diag- 
onals; 

(2) V is invertible; 

(3) N is nilpotent: If d - 1 = dim (M^) then N"^-^ = 0; 

(4) [A^, D] = ND - DN = 0; 

(5) G = V{D + N) 

Let X e M^, and let n > d. Using (2)-(5), we get a-"G"x = Va" {D + A^)" V~^x 
YliZo {^)Va~'^D'^~^N^v^^x. But the matrix a^"D"~* is diagonal with entries 
{a~^s^~'^\s G spec (F)\{a}} , < i < d — 1. Using finally assumption (ii), 
we conclude that 



lim ( ^ ) a""s"-^ = 0, (69) 

n— >c« \i J 

and the proof of fl67|) is completed. ■ 

Proof. (Remark II V.3|) Let the conditions be as stated in the Remark. From 
the arguments in the proof of Lemma IIV.21 we see that the two vectors on 
the left-hand side in ( l60l) may be decomposed as follows: 

a-"F"x = {w\x)w + (1 - a""G") (a - G)'^ t] + a~''G''PM^x (70) 

and 

{'w\x)^ = {'w\x)'w + {a — G)~^ T]. (71) 
Hence, the difference is in A^-*-, and 

||a-"F"a; - {w \x\)i\\ = \\a-'"G" {Pm^x - (a - G)'^ r]) \\ 

<Cn maxj- |s G spec (F)\{a}| 

which is the desired conclusion. (1601) . ■ 
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